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ABSTRACT 

We present optical spectroscopy obtained with the Space Telescope Imaging Spec- 
trograph (STIS) of five young massive star clusters in the starburst galaxy M82. A 
detailed analysis is performed for one cluster 'M82-A1' and its immediate environment 
in the starburst core. From HST archive images, we find that it is elliptical with an 
effective radius of 3.0 ± 0.5 pc and is surrounded by a compact (r = 4.5 ± 0.5 pc) 
H II region. We determine the age and reddening of M82-A1 using synthetic spectra 
from population synthesis models by fitting both the continuum energy distribution 
and the depth of the Balmer jump. We find an age of 6.4 ± 0.5 Myr and a photo- 
metric mass estimate of M = 7-13 x 10 5 M©. We associate its formation with the 
most recent starburst event 4-6 Myr ago. We find that the oxygen abundance of the 
Hll region surrounding M82-A1 is solar or slightly higher. The Hn region has a high 
pressure P/k = 1-2 xlO 7 cm _3 K. The diffuse gas in region A has a slightly lower 
pressure, which together with the broad Ha emission line width, suggests that both 
the thermal and turbulent pressures in the M82 starburst core are unusually high. 
We discuss how this environment has affected the evolution of the cluster wind for 
M82-A1. We find that the high pressure may have caused the pressure-driven bubble 
to stall. We also obtain spectroscopic ages for clusters Bl-2 and B2-1 in the 'fossil' 
starburst region and for the intermediate age clusters F and L. These are consistent 
with earlier studies and demonstrate that star formation activity, sufficiently intense 
to produce super star clusters, has been going on in M82 during the past Gyr, perhaps 
in discrete and localized episodes. 

Key words: galaxies: evolution - galaxies: individual: M82 - galaxies: ISM - galaxies: 
starburst - galaxies: star clusters. 



1 INTRODUCTION 

Young, massive super star clusters (SSCs) 1 are a prominent 
product of the intense star formation episodes that occur in 
nearby starbursts, and presumably in merging galaxies at 
high redshift. They are important objects to study because 
of the insights they provide on: the formation and destruc- 
tion of globular clusters; star formation processes in extreme 
environments; and the triggering and feeding of supergalac- 
tic winds. The properties of young massive star clusters are 



1 In this paper, we define SSCs to include star clusters with half- 
mass radii of ^ 5 pc; ages of 500 Myr; and total masses of 
^ 10 5 M 



extensively reviewed in a recent proceedings (Lamers, Smith 
& Nota 2004). 

M82 is the best nearby example of a giant starburst 
galaxy at a distance of 3.6 Mpc (Freedman et al. 1994). 
The present intense burst of nuclear star formation was trig- 
gered by a tidal interaction with M81 some 2-5 x 10 8 yr ago 
(Brouillet et al. 1991; Yun, Ho & Lo 1994). The active star- 
burst region has a diameter of 500 pc and is defined optically 
by the high surface brightness regions or clumps denoted 
A, C and E by O'Connell & Mangano (1978). An optical 
imaging study of the central region of M82 with the Plan- 
etary Camera onboard the Hubble Space Telescope (HST) 
by O'Connell et al. (1995) revealed the presence of over one 
hundred candidate SSCs within the visible starburst. They 
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find that the optically-bright structures containing the clus- 
ters correspond to known sources at X-ray, infra-red (IR) 
and radio wavelengths. They argue that this cannot be ac- 
cidental and suggest that regions A, C and E represent the 
parts of the starburst core which are the least obscured along 
the line of sight. 

Region A is of special interest. It contains a remark- 
able complex of SSCs with very high continuum and emis- 
sion line surface brightnesses (O'Connell & Mangano 1978; 
O'Connell et al. 1995). The well-known large-scale bipolar 
outflow or 'superwind' appears to be centred on regions A 
and C (Shopbell & Bland-Hawthorn 1998; Gallagher et al., 
in prep.). In a recent paper, Melo et al. (2005) have cata- 
logued the young stellar clusters in the starburst core of M82 
using HST archive data. They find a total of 197 optically 
visible clusters with 86 of these residing in region A. 

The age and evolution of the starburst in M82 have 
been the subject of numerous studies. The starburst region 
appears to be composed of two spatially separate stellar pop- 
ulations. The presence of a strong IR continuum and large 
CO absorption index indicates a population of red super- 
giants in the core (Rieke et al. 1993; Satyapal et al. 1997; 
Forster Schreiber et al. 2001). The youngest and most mas- 
sive stars responsible for the strong Lyman continuum flux 
are located in a ~ 200 pc circumnuclear arc (Achtermann & 
Lacy 1995). Using evolutionary synthesis models, Rieke et 
al. (1993) and Forster Schreiber et al. (2003) find that two 
bursts are needed to reproduce the observational properties. 
The latter authors find the bursts occurred « 10 and 5 Myr 
ago with a duration of a few million years. The first burst 
was particularly intense in the core and the second occurred 
mainly in the circumnuclear regions. 

While models of the active starburst region suggest star 
formation is propagating outwards, studies of star clusters 
outside the central regions reveal considerably older star for- 
mation episodes. Gallagher & Smith (1999) derived ages of 
60±20 Myr for two SSCs M82-F and L located 440 pc south- 
west of the nucleus. Photometric age-dating of the extended 
"fossil" starburst region B, located 0.5-1 kpc north-east 
of the nucleus, shows that the peak of the star formation 
episode occurred at log(age/yr) = 9.0 ± 0.4 (de Grijs et al. 
2001, 2003, 2005). 

One long outstanding question concerning starburst en- 
vironments is whether the initial mass function is abnormal 
in the sense that less low mass stars are produced in an in- 
tense star-forming episode. M82 has been the focus of many 
studies aimed at settling this question. The mass function of 
the central starburst region of M82 has been estimated based 
on the light-to-mass ratio, but this is controversial mainly 
owing to the difficulty of determining extinction corrections 
in such a complex system. Some studies (e.g. O'Connell & 
Mangano 1978; Telesco et al. 1991; Satyapal et al. 1997) 
find a normal mass function, while others (e.g. Rieke et al. 
1993; Forster-Schreiber et al. 2003) find a deficiency of low 
mass stars. Light-to-mass ratios determined for individual 
SSCs are more robust. Smith & Gallagher (2001) measured 
the dynamical mass of M82-F and showed that its visual 
mass-to-light ratio required a mass function deficient in low 
mass stars, and thus it is unlikely to become an old globu- 
lar cluster. McCrady, Gilbert & Graham (2003) determined 
dynamical masses for two SSCs on the western edge of the 
active starburst region in M82, and found that one object 



appeared to be deficient in low mass stars. A re-analysis 
of M82-F by McCrady, Graham & Vacca (2005) based on 
HST images obtained with the Advanced Camera for Sur- 
veys (ACS) and Keck/NIRSPEC spectra confirm that it has 
an abnormal mass function, and thus possibly a peculiar 
stellar mass function. 

In contrast, analyses of nebular line diagnostic ratios 
indicate a lack of hard-UV photons in some starburst re- 
gions, implying a deficit of high mass stars (e.g. Puxley et al. 
1989; Achtermann & Lacy 1995). In a recent study, Rigby & 
Rieke (2004) suggest an alternative explanation for the low 
[Ne lll]/[Ne n] ratios observed in solar metallicity-type star- 
bursts. They propose that the hard ionizing photons from 
massive stars are absorbed in long-lived ultracompact Hn 
regions which are a consequence of the high interstellar den- 
sities encountered in a starburst. 

The proximity of M82 (0".l = 1.8 pc) means that in- 
dividual clusters are easily resolved by HST and are bright 
enough to obtain spectroscopy. In this paper, we present 
spectroscopy obtained with the Space Telescope Imaging 
Spectrograph (STIS) onboard HST for five massive clus- 
ters of varying ages. The first object, which we designate 
as M82-A1, is a bright isolated cluster in region A. We will 
focus on the analysis and interpretation of this STIS dataset 
because it represents the first optical spectroscopy of an in- 
dividual cluster and its immediate environment in the star- 
burst core. We also present STIS spectroscopy for one slit 
position crossing both clusters F and L, and slits covering 
the two brightest clusters in the fossil region B. These spec- 
tra are of lower quality than those of M82-A1 but we are 
able to derive spectroscopic ages. 

The STIS slit position for cluster M82-A1 extends 
across regions A and C, close to the nucleus and the base 
of the superwind, sampling the ionized gas in the starburst 
core on spatial scales of a few parsecs. In a separate paper 
(Smith et al., in prep.), we will discuss the nature of the ion- 
ized component in the M82 starburst, as revealed by these 
optical emission line data. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Spectra 

We obtained long-slit HST/STIS spectra at four positions 
in the centre of M82 (GO 9117; P.I. O'Connell). One slit 
was centred on the brightest isolated cluster in region A, 
which we designate as M82-A1. Another slit sampled M82- 
F and M82-L, while the two remaining pointings were cen- 
tred on M82-B1-1 and M82-B2-1, within the post-starburst 
region (de Grijs et al. 2001). In Fig. we show the slit 
positions for the five clusters superimposed on an archive 
HST+ACS/WFC F814W image (see Sect.l2~2l. 

Details of the STIS spectroscopy observations are listed 
in Table Spectra were taken in three different bands: the 
blue continuum (including the Balmer jump) , the Ha region, 
and the near-infrared continuum. Coordinates of the four 
clusters initially placed at the centre of the slit are given in 
the footnote to the table. In each case, the position angle was 
oriented to place other interesting regions on the slit. In the 
case of Al and F+L, multiple exposures were made, stepped 
along the slit to improve the elimination of hot pixels in the 
data reduction process. 
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Figure 1. HST+ACS/WFC F814W image of M82 showing the slit locations for the star clusters Al, F, L, Bl-1 and B2-1. The more 
prominent cluster complexes (A, B, C and E) arc labelled. The inset in the lower left-hand corner (corresponding to the white outline) 
shows the position of the STIS slit across region A and the location of cluster Al. 
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Table 1. HST/STIS Spectroscopic Observations 



Grating 



Cluster 


bllt 


Name 

A range (A) : 
AA (A pix J : 


G430L 

2900-5700 

2.73 


G750M 

6295-6865 

0.56 


G750L 

5250-10300 

4.92 


Al 


Ap=52" x O'.'l 


No. of exposures: 


12 


4 


— 




PA= 229° 


Date : 


2002/11/27 


2002/11/26 






Oflfsets=Multiple 


Exposure time : 


231.4m 


70.8m 




Bl-1 


Ap= 52" x O'.'l 


No. of exposures: 


3 


2 






PA= 46° 


Date : 


2002/06/03 


2002/06/03 






Offsets=none 


Exposure time : 


52.2m 


39.m 




B2-1 


Ap= 52" x O'.'l 


No. of exposures: 


3 


2 






PA= 235° 


Date : 


2002/11/02 


2002/11/02 






Offsets=none 


Exposure time : 


52.2m 


39.1m 




F+L 


Ap= 52" X O'.'l 


No. of exposures: 






4 




PA= 207° 


Date : 






2001/12/26 




Offsets=multiplc 


Exposure time : 






163.7m 



Al : a = 09 h 55 m 53 s . 42; <5 = +69° 40' 51". 08 (J2000) 
F : a = 09 h 55 m 47 s . 01; <5 = +69° 40' 42". 30 (J2000) 
Bl-1: a = 09 h 56 m 03 s . 34; <5 = +69° 41' 12". 20 (J2000) 
B2-1: a = 09 h 55 m 54 s . 55; 5 = +69° 41' 01". 60 (J2000) 
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Figure 2. STIS G750M two-dimensional spectral image showing the cluster M82-A1 and its compact Ho region over the wavelength 
range 6500-6800A. The y-scalc is in arc seconds and the slit orientation is shown. The nebular emission lines are marked and the surface 
brightness scale is in units of 10~ 16 erg cm~ 2 s _1 A -1 arcsec~ 2 . 
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Figure 3. G430L spectra of the clusters Bl-1 and B2-1 (left-hand side) and G750L spectra of the clusters F and L (right-hand side). 
Each spectrum has been smoothed by 2 A and the main stellar absorption features are labelled. The flux is in units of 10 — 16 erg cm -2 



To combine the individual image sets of Al, Bl-1 and 
B2-1, we first registered the cosmic ray-corrected images 
from the calstis data reduction pipeline. We then flipped 
the cold pixels present in the datasets to high values by in- 
specting intensity histograms of the images. We next used 
the iraf/stsdas cosmic ray rejection (occrej) task to com- 
bine all the images and eliminate the remaining cosmic rays 
and hot and cold pixels. Finally we used the x2d task to pro- 
vide rectified absolute flux and wavelength calibrated two- 
dimensional spectral images. 

In Fig. |5J a small portion of the reduced two- 
dimensional image of M82-A1 for the G750M grating is 
shown. The cluster is clearly resolved and is surrounded by a 
compact H n region. From cuts across the STIS 2D spectra, 
we measure a FWHM of ~ 6 pixels or 0".3 for M82-A1, ex- 
ceeding the slit width of 0".l. The spectral resolution of the 
G430L and G750M gratings for an extended source is 2-3 
pixels. We measure a resolution of 6. 7 ±0.6 A from Gaussian 
fits to the unresolved nebular H/3 emission line originating 
in ionized gas along the slit. We will therefore assume for 
this cluster that the actual spectral resolution is 2.5 pixels 
or 6.7 A (G430L) and 1.4 A (G750M grating). 

To extract one-dimensional spectra, we used an extrac- 
tion width of 14 pixels for clusters Al, Bl-1 and B2-1, to 
include all the cluster light. For Al, we performed a back- 
ground subtraction using 10 pixel-wide regions centred at 
±25 pixels from the centre of the cluster. These particular 
regions were chosen because they contain negligible emis- 



sion from the ionized gas in region A compared to the clus- 
ter (< 5% at Ha; cf. Fig|5J. The exact choice of background 
bins for the clusters in Bl and B2 was not critical, so we used 
a prescription analogous to that for region A. No correction 
for slit losses have been made. 

Clusters F±L were observed with the G750L grating. 
To enable the correction of fringing effects, which are par- 
ticularly prominent for redder wavelengths, the science ex- 
posures were interspersed with CCD flatfields. The obser- 
vations were calibrated using the CALSTIS pipeline within 
iraf/stsdas. After basic reduction and cosmic ray rejec- 
tion, normalized fringe flats were prepared for each exposure, 
following the procedure outlined in STIS/ISR98-29. These 
fringe flats were then applied to correct the fringing in the 
spectra. Careful examination of the fringe-corrected spectra 
showed this procedure to be quite clean and accurate. In 
this way, we reduced the dataset to four, cosmic-ray cleaned 
and defringed, two-dimensional spectra. 

Before the process of final calibration, during which 
a spectrum is rectified along the wavelength and spatial 
axes, as well as converted to physical flux units, the spec- 
tra from each of the four exposures needed to be registered, 
shifted and combined into a single two-dimensional spec- 
trum, with the associated removal of hot pixels. The com- 
bination was performed as follows: first, the spatial profile 
of the cluster was obtained by summing over wavelength 
pixels in bins which were free of contamination from both 
emission lines and hot pixels lying near the position of the 



6 L.J. Smith et al. 



Table 2. HST images and derived parameters for cluster M82-A1. The errors in the effective radius and 
magnitude measurements are ±30 mas and ±0.10 mag respectively. 



Filter 


Camera 


Plate Scale 


Exposure 


FWHM 


Minor/Major 


_R c ff 


VEGAMAG 




± CCD 


(" Pix- 1 ) 


Time 


(major axis) 


Axis Ratio 




Photometry 








(s) 


(mas) 




(mas) 


( ma.er^ 


F439W 


WFPC2/WF4 


0.099 


3300 


140 


0.74 


130 


18.77 


F439W 


WFPC2/WF4 


0.099 


4100 


140 


0.71 


110 


18.75 


F555W 


WFPC2/WF4 


0.099 


2500 


170 


0.68 


160 


17.23 


F555W 


WFPC2/WF4 


0.099 


3100 


170 


0.65 


160 


17.22 


F547M 


WFPC2/PC 


0.045 


100 


180 


0.72 


180 




F814W 


WFPC2/WF4 


0.099 


140 


220 


0.66 


210 


14.77 


F814W 


ACS/WFC 


0.049 


120 


210 


0.70 


210 


14.91 


F160W 


NICM0S/NIC1 


0.043 


448 


200 


0.79 


200 




F190N 


NICMOS/NIC2 


0.076 


384 


200 


0.71 


190 




F222M 


NICMOS/NIC2 


0.076 


448 


160 


0.79 


160 




F656N 


WFPC2/PC 


0.045 


1000 


290 


0.57 


260 


8.54E-14 1 



1 Total Ho flux in units of ergs cm 2 s 1 



cluster in the spectrum. A separate profile was made for 
each exposure. The peak of the cluster was determined by 
fitting gaussians to each of the four cluster profiles. The 
two-dimensional spectra were then shifted along the spatial 
axis by the amount needed to place the peak of the cluster 
exactly at the central spatial pixel of the spectrum. These 
shifted spectra were finally combined using the STSDAS rou- 
tine mscombine, which cleanly removes hot pixels from the 
final spectrum. 

Since the clusters are marginally resolved with the 0.1" 
aperture of the slit, aperture illumination effects can be sub- 
stantial. Therefore, we used the STSDAS routine xlD, part of 
the CALSTIS ensemble, to extract the final one-dimensional 
spectrum of clusters F and L within 11 pixel wide bins, as 
well as adjacent galaxy backgrounds. Using XlD ensures 
that wavelength dependent aperture illumination effects are 
taken into account during flux calibration of the spectrum. 
The backgrounds were then subtracted from the cluster ex- 
tractions to obtain spectra of pure cluster light. For cluster 
F, the background accounts for no more that 10% of the light 
in the cluster extraction, but this fraction is higher for the 
lower luminosity cluster L, especially in the blue end of the 
spectrum, since cluster L is heavily reddened by intervening 
dust. 

In Fig. [3]we show the final reduced G430L spectra for 
clusters Bl-1 and B2-1, and the G750L spectra for clusters F 
and L. The spectra for cluster Al are presented in Sect. 13.21 



2.2 Images 

In this section we present HST archival images for cluster 
M82-A1 and measure its size and perform photometry. HST 
images and measurements have been presented by de Grijs 
et al. (2001) for clusters Bl-1 and B2-1. In their nomen- 
clature, Bl-1 corresponds to No. 28 in their table 2, and 
B2-1 to No. 12 in their table 3. The V magnitudes of Bl- 

1 and B2-1 are given as 17.93 and 17.90 mag respectively. 
Smith & Gallagher (2001) present HST images of clusters F 
and L obtained with the Wide Field and Planetary Camera 

2 (WFPC2) while McCrady et al. (2005) present images of 



cluster F obtained with the High Resolution Channel (HRC) 
of the ACS instrument. 

We obtained HST images taken with the WFPC2, ACS 
and NICMOS cameras for the central region of M82 from the 
HST archive; the details of these images are given in Tabled 
The broad-band WFPC2 F439W, F555W and F814W im- 
ages from the GO 7446 programme (PI O'Connell) are cen- 
tred on the fossil starburst regions Bl and B2 (see de Grijs 
et al. 2001) and region A is on the WF4 CCD. We also ob- 
tained narrow-band images in the F656N and F547M filters 
which are part of the GO 6826 programme (PI Shopbell) on 
the superwind of M82. We used the F547M filter for con- 
tinuum subtraction in the F656N filter. These images are 
discussed in more detail in Gallagher et al. (in prep.). The 
single ACS F814W image was obtained as part of the GO 
9788 snapshot programme (PI Ho) and is shown in Fig. 
The NICMOS images from the GO 7218 programme (PI 
Rieke) are described in Alonso-Herrero et al. (2003). 

Cosmic rays were removed from the pipeline-reduced 
images using the Laplacian cosmic ray identification routine 
of van Dokkum (2001). The individual images were then 
checked for saturation and averaged. Photometry and radius 
measurements of M82-A1 are complicated by the varying 
background of emission and absorption by gas and dust in 
the starburst core. The images also show that the cluster is 
quite elliptical (Fig. . 

To measure the size of M82-A1, we used the ishape 
algorithm (Larsen 1999) together with the tinytim pack- 
age (Krist 2004) to correct for the point-spread function 
(PSF) for the WFPC2 and NICMOS cameras. For the ACS 
F814W image, we used an empirical PSF generated from 
point sources in the field of 47 Tuc (Mora, priv. comm.). 
First, we determined the best fitting radius for the profile 
fits by examining fits to the F547M image using a Moffat 
function with a power index of 1.5. We found that a radius 
of 0.5 arcsec was optimal since this was sufficiently large 
but excluded light from other nearby sources. We then ex- 
perimented with different profile types and chose a Moffat 
function with a power index of 1.5 since this gave the lowest 
residuals in the overall fit. We also tried varying the power- 
law index of the Moffat profile but found that a value of 1.5 
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consistently gave the best fit. Elson, Fall & Freeman (1987) 
determined that profiles of this type provide good fits to 
LMC star clusters while de Grijs et al. (2001) found similar 
results for the profiles of the brightest clusters in region B 
of M82, including Bl-1 and B2-1. 

In Table|5J we list the FWHM of the major axis and the 
minor/major axis ratio, as estimated by ishape for each im- 
age. We also give the effective or half-light radius R c g where 
we have used the expressions given by Larsen (2004) to con- 
vert the FWHM of an elliptical Moffat profile fit to a half- 
light radius. We find very consistent results for the WF4 and 
ACS F814W measurements given the different pixel scales 
and PSFs. For comparison, as detailed in the instrument 
handbooks, the PSFs of the WF4, PC and ACS/WFC chips 
are w 80, 130 and 80 mas respectively. To estimate appro- 
priate errors, we fitted the F547M and F814W/ACS im- 
ages with a number of different profile types. We find that, 
in terms of the residuals, King profiles with different con- 
centration parameters also give acceptable fits. From the 
scatter in the resulting R c g values, we find that the likely 
error in each measurement is ±30 mas. Taking a straightfor- 
ward mean of the radii measured in the continuum filters, 
we find i? cfi = 170 ± 30 mas or 3.0 ± 0.5 pc for M82-A1. We 
derive a mean minor/major axis ratio of 0.72 ± 0.05 from 
the same images. For comparison, McCrady et al. (2003) 
have measured R e a for several clusters (including M82-A1 
which is MGG-al in their designation) from NICMOS NIC2 
F160W images. From spherical King profile fits, they find 
R c h = 119± 14 mas which is smaller than our own more pre- 
cise measurements which allow for ellipticity. It is apparent 
from Table [5] that the radius of M82-A1 decreases towards 
shorter wavelengths; we discuss this effect in Sect. 13. ll 

Measurements of the size and total flux of the com- 
pact Hn region surrounding M82-A1, based on Ha emission 
measured in the F656N filter, are also given in Table [5] It is 
highly elliptical with a major axis FWHM value of 290 mas 
and an axis ratio of 0.57, giving a half-light radius R e ff of 
260 ±30 mas or 4.5 ±0.5 pc. Thus the Hn region is only 1.5 
times larger than the cluster. 

In Table we list the magnitudes of M82-A1 for the 
wide-band filters in the VEGAMAG system. M82-A1 is su- 
perimposed on a varying background. We tried various ap- 
proaches to subtracting the background and found that the 
most robust and simplest method was to use a flat plane. 
The level was calculated from an average of measurements, 
taken at four locations around the cluster M82-A1, that 
were deemed to contain "representative" background val- 
ues. The photometry of M82-A1 was performed using the 
PHOT task within the NOAO digiphot package. It was de- 
cided that a 7 pixel radius aperture gave the best results, 
considering the morphology of M82-A1 is such that it has 
a protrusion extending north-westwards of the cluster. The 
total counts measured in each aperture were corrected for 
CCD charge transfer effects using the formulation of Dol- 
phin (2002). They were then converted to magnitudes us- 
ing the zeropoints provided for the VEGAMAG system in the 
instrument handbooks. We estimate the measurements are 
accurate to ±0.10 mag because of the uncertainties in the 
varying background. This level of accuracy was also verified 
with two sets of independent measurements. 
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Figure 4. Measured effective radii (in milli-arcsec) of M82-A1 
plotted as a function of wavelength using the data in Table [5] 
The errors bars represent the measurement uncertainties of ±30 
mas. 



3 CLUSTER M82-A1 

3.1 The change of radius with wavelength 

Before presenting and discussing the STIS spectra of M82- 
Al, we first discuss possible reasons for the decrease in ra- 
dius of the cluster at progressively shorter wavelengths. In 
Fig.0]we plot R c s against wavelength for M82-A1. If we ex- 
clude the NICMOS measurements in Table [5] (because the 
diffraction limit is close to the radius of M82-A1), then the 
cluster radius appears to decrease at successively shorter 
wavelengths; the radius is 30% smaller in the F439W filter 
compared to the F814W filter. 

McCrady et al. (2005) find that the radius of cluster 
M82-F decreases with increasing wavelength, and they in- 
terpret this as being due to mass segregation, with the more 
centrally concentrated stars being the most massive cooler 
evolved stars. For M82-A1, we see the opposite effect with 
the radius deceasing towards blue wavelengths. 

It is not clear what is causing this decrease in radius. 
At the young age of this cluster, most of the detected light 
should originate from stars with roughly the same mass. We 
find from the output of the spectral synthesis code Star- 
burst99 (see Sect. 13. 31 that this is indeed the case, and that 
any change in colour with mass will be shallow for a mass 
segregated cluster. This leads us to conclude that it is diffi- 
cult to produce a radial gradient due to mass segregation. 

We now consider whether the compact Hn region sur- 
rounding the cluster could affect the radius measurements. 
It is unlikely that contamination from nebular emission lines 
in the broad-band filters is causing the radius variations be- 
cause the emission lines are weak (see Figs.|K|and |SJ and the 
F814W filter does not contain any significant nebular emis- 
sion. The presence of dust in the Hn region may, however, 
produce a decrease in radius towards bluer wavelengths. A 
diffuse foreground screen of dust will change the brightness 
but, since it will dim the cluster uniformly, it will not affect 
the radius or shape as a function of wavelength. However, if 
dust is concentrated in the H II region, then the cluster would 
maintain a roughly constant shape (as measured), yet be 
dimmer and smaller with decreasing wavelength because of 
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Table 3. Absorption and emission line measurements for M82-A1. Velocities are in the heliocentric frame 
of reference. 



Line 


■^vac 


Velocity 


FWHM 


Equivalent 


Comments 




(A) 


(km s -1 ) 


(km s -1 ) 


Width (A) 




G430L spectrum 










[On] 


3728.48 


448 




-2.9 ± 0.6 


blended doublet 


H8 abs 


3890 1 5 


"1 "71 
1 / 1 




3.2 ± 0.5 


blended with He 1? 


Can K 


3934.78 


205 




3.8 ± 0.6 


Interstellar 


Can H 


3969.59 


415 






Interstellar; blended with He 


He abs 


3971.19 


294 




3.9 ± 0.7 


blended with Call H 


H5 abs 


4102.89 


335 




3.0 ±0.7 




H7 em 


4341.68 


332 








H7 abs 


4341.68 


327 




2.0 ±0.9 




H/3 em 


4862.68 


298 




-1.8 ±0.2 




[Om] 


4960.30 


343 




-1.0 ± 0.1 




[Oiii] 


5008.24 


319 




-1.8 ± 0.2 




G750M 


spectrum 










[Nn] 


6549.85 


322 


9.3 


-4.8 ± 0.2 




Ha 


6564.61 


319 


95 


-27.6 ± 0.3 




[Nn] 


6585.28 


318 


96 


-15.5 ± 0.2 




[S 11] 


6718.29 


319 


77 


-2.3 ± 0.2 




[Sn] 


6732.67 


322 


84 


-3.1 ± 0.2 





the greater obscuration at shorter wavelengths. We there- 
fore conclude that the most probable explanation for the 
decrease in radius with decreasing wavelength is the effect 
of dust in the surrounding compact Hn region. To resolve 
this conclusively, ACS/HRC observations of M82-A1 are re- 
quired. Dynamical models addressing the colour gradients in 
young mass segregated clusters and a check for cool super- 
giants that could support colour gradients are also needed. 



3.2 Spectra 

In Figs.|S]and|S] we show the flux-calibrated spectra of M82- 
Al for the G430L and G750M gratings. The major line iden- 
tifications, heliocentric velocities, FWHMs (for the G750M 
spectrum only) and equivalent widths are presented in Ta- 
ble [3 for the absorption and emission features. The equiva- 
lent widths were measured by defining continuum windows 
in the spectra and the errors represent la uncertainties 
based on the error in the continuum. 



3.2.2 Emission Lines 

A number of emission lines arising from the compact Hn 
region surrounding M82-A1 are apparent in the STIS spec- 
tra. We detect nebular emission lines of Ha, H/3, H7, [O n] 
A3727,3730, [Om] A4960, 5008, [Nn] A6550,6585 and [Sn] 
A6718, 6733. The low [S ll]/Ha ratio of 0.21 is consistent with 
photoionization rather than shock heating. In Fig. |S| the 
red member of the [S 11] doublet is the strongest of the two 
lines, indicating that the Hn region has a high density. In 
Section r3.4l we derive the electron density and oxygen abun- 
dance. 

The comparatively high resolution of the G750M spec- 
trum allows us to measure accurate velocities and FWHMs 
for the H 11 region, as given in Table 2. We find a mean he- 
liocentric velocity for the Ha, [N n] and [S n] nebular lines 
of 320 ±2 krns" 1 , in excellent agreement with the measured 
value of the cluster radial velocity. The lines have an ob- 
served FWHM of 89 ± 8 kms" 1 and are clearly resolved 
since the measured resolution is 63 kms -1 ; the intrinsic 
mean FWHM is then 62 ± 11 kms - . This gives a maxi- 
mum expansion velocity for the Hn region of ~ 30 kms -1 . 



3.2.1 Absorption Lines 

The G430L spectrum is dominated by Balmer absorption 
lines from the cluster, and a weak Balmer jump at 3650 A, 
indicating a young age. Strong interstellar Call absorption 
lines are also present. No stellar absorption lines are detected 
in the G750M spectrum. The mean heliocentric radial veloc- 
ity of the cluster as measured from the Balmer absorption 
lines is 319 ±22 kms - . The velocity of the interstellar Can 
K line is ±205 kms -1 . This is very close to the systemic 
velocity of M82 of ±200 kms -1 (Achtermann & Lacy 1995). 
The cold interstellar gas therefore appears to be at rest for 
this line of sight near the nucleus of M82. 



3.3 Cluster parameters 

The starburst core of M82 suffers from a high and non- 
uniform extinction. For example, Satyapal et al. (1995) 
find values for the visual extinction Av ranging from 2- 
12 mag from near-IR nebular line ratios, assuming a fore- 
ground screen model. For the whole M82-A complex, centred 
about 6" SW of M82-A1 (see Figs. [T]i and[TJ>), Marcum & 
O'Connell (1996) estimated that E(B - V) ~ 0.60 arises in 
a foreground dust screen with an additional r(5500 A) ~ 0.5 
arising from internal dust. As can be seen from Figs. Q and 
[5] the STIS slit crosses some areas of complete obscuration 
at optical wavelengths. The cluster M82-A1 sits close to the 
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Figure 5. G430L spectrum of the cluster M82-A1 and its Hn region smoothed by 2 A. The nebular emission and cluster absorptii 
lines are marked above and below the spectrum respectively. The flux is in units of 10 — 15 erg s _1 cm -2 A -1 
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Figure 6. STIS G750M spectrum of the cluster M82-A1 and its Hn region. The nebular emission lines from the surrounding Hn region 

-2 A -l 



conspicuous dust lane that bisects the galactic disk; to the 
east of this cluster (or lower pixel number in Fig. 0, we 
detect little emission due to the very high extinction. 

Ages and reddening for clusters less than 10 Myr old are 
usually determined from spectroscopy using the UV wave- 
length region because the spectral energy distribution (SED) 
is largely independent of age and the stellar wind features 
are sensitive age indicators (e.g. Tremonti et al. 2001; Chan- 
dar et al. 2004). However, given sufficient spectral resolution, 



the short-wavelength optical band (especially the continuum 
structure and strong absorption lines in the Balmer jump re- 
gion between 3500 and 4100 A) should provide a comparable 
hold on age and reddening. Gallagher & Smith (1999) de- 
rived the age of cluster M82-F by comparing the observed 
spectrum in the region of the Balmer jump with theoretical 
model cluster spectra. Here, we use recent high-resolution 
population synthesis models to fit the combined spectrum 
of M82-A1. 
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Figure 7. Comparison of the de-reddened (E(B - V = 1.28)) STIS spectrum of M82-A1 with a synthetic Starburst99 (SB99) spectrum 
for an age of 6.4 Myr. For this combination of age and reddening, the depth of the Balmer jump and the overall SED are well-fitted. The 




3600 4000 4400 

Wavelength (A) 

Figure 8. Comparison of the de-reddened STIS spectrum of M82- 
Al with synthetic SB99 spectra for ages of 5.5 Myr (E(B — V) = 
1.42) and 7.5 Myr (E(B - V) = 0.95) demonstrating that the 
depth of the Balmer jump cannot be fitted at these ages. The 
flux is in units of 10 — 14 erg s — 1 cm -2 A -1 and the spectra for 
the 7.5 Myr comparison have been multiplied by a factor of 2.5. 



First, we binned the G750M STIS spectrum to match 
the G430L spectrum and merged the two spectra. We fit 
the combined spectrum using v5.0 of the spectral synthe- 
sis code Starburst99 (SB99: Leitherer et al. 1999; Gonzalez 
Delgado et al. 2005) which includes a new synthetic high res- 
olution (0.3 A) spectral library from 3000-7000 A containing 
non-LTE models for OB stars (Martins et al. 2005). We gen- 
erated a series of SB99 models for different ages (2-10 Myr), 
assuming solar metallicity (see Sect. 13.41 . an instantaneous 
burst with a Kroupa (2001) IMF, lower and upper masses of 
0.1 and 100 M©, and the enhanced mass loss Geneva tracks. 



The resulting high resolution synthetic spectra were 
smoothed, binned, and velocity shifted to match the G430L 
spectrum of M82-A1. We de-reddened the observed com- 
bined spectrum using a foreground dust screen, employ- 
ing the Galactic reddening law from Howarth (1983) and 
R = A(V)/E(B~V) =3.1, and a range of E(B~V) values 
from typically 0.9-1.8 mag. These de-reddened spectra were 
then normalised to the synthetic spectra at 6800 A and the 
quality of the fit for each possible combination of age and 
E(B — V) was judged by the x 2 statistic, with highest weight 
given to the spectral region 3700-4400 A. Typically, we find 
that we can match the model continuum longward of the 
Balmer jump from 3700-6800 A by varying E(B - V) until 
a good fit is obtained; this then provides a value of the red- 
dening for each model as a function of age. The quality of 
the fit of the model spectrum to the observed de-reddened 
spectrum below the Balmer jump then determines the age. 

Using this approach, we find that the STIS spectra 
of M82-A1 can only be matched by the SB99 synthetic 
spectra for a relatively narrow age range of 5-7 Myr and 
E(B~V) = 1.5-1.2 mag. We checked the SB99 results with a 
comparable grid taken from the independent synthesis mod- 
els of Bruzual & Chariot (2003; BC03), which incorporate 
a lower resolution (3 A) but more realistic empirical spec- 
tral library. The fits are quite similar to those of SB99, with 
a slightly larger age range of 4-7 Myr. From ground-based 
spectrophotometry, Marcum & O'Connell (1996) obtained 
an average age of ~ 5 Myr for the entire M82-A complex, in 
good agreement with these estimates for M82-A1. 

In Fig. □ we show the SB99 fit for the best fitting SB99 
model, at an age of 6.4 Myr and E(B - V) = 1.28 mag. The 
overall fit to the continuum is good, and it can be seen that 
the upper members of the Balmer series are also fairly well- 
fitted (e.g. H8 which has minimal nebular line emission). 
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This model, which has the largest Balmer jump in either 
model grid, slightly overestimates the observed flux in the 
Balmer continuum below 3600 A. The discrepancy could be 
caused by internal extinction within M82-A1 or perhaps a 
small error in the STIS calibration here. Alternatively, there 
could be weaknesses in the models. Adjacent models which 
fit the continuum above the Balmer jump by varying the 
E(B — V) parameter provide poorer short-wavelength fits 
below the Balmer jump, as illustrated in Fig. |S] Abrupt 
changes in the Balmer jump and the color of the contin- 
uum in the SB99 models for ages of 6-10 Myr are probably 
caused by intermediate temperature class I-II stars, and it is 
not clear that these are being accurately modelled (Gonzalez 
Delgado et al. 2005). 

We can now employ two independent indicators for red- 
dening and age to determine if they are consistent with 
those found from the spectral analysis. First, we use the 
nebular H I emission lines to measure the interstellar extinc- 
tion by comparing the observed Hcv/H/3 and H7/H/3 ratios 
with the intrinsic case B values of Hummer & Storey (1987). 
We correct for stellar absorption by directly subtracting the 
reddened synthetic SB99 spectrum from the observed STIS 
spectrum using the reddening values we have determined for 
the best-fit age range. We find E{B - V) = 1.38 ± 0.15 mag, 
in good agreement with dereddening the cluster energy dis- 
tribution. This method of directly subtracting the synthetic 
spectrum to recover the intrinsic nebular line fluxes has the 
advantage that we can recover the H7 emission line and ob- 
tain a more accurate measurement for the H/3 line flux. We 
also compare the observed (B — V) and (V — I) colours of 
M82-A1 with the predicted colours from SB99 for ages of 6 
and 7 Myr and find E(B - V) = 1.54 ± 0.26 in reasonable 
agreement with the above. For the reddening towards M82- 
Al, we will adopt the mean value given by the cluster SED 
and the nebular lines of E(B — V) = 1.35 ± 0.15 mag. 

Our independent age indicator is the strength of the ob- 
served equivalent width of the nebular Ha emission (Table[3Jl 
when compared with the predicted values from Starburst99. 
The observed equivalent width of 32 ± 2 A (corrected for 
stellar absorption) gives an upper limit to the age of 6.7 
Myr. This value is consistent with the age determined from 
the Balmer jump, and suggests that few ionizing photons 
escape from the Hn region. We thus conclude that the age 
of M82-A1 is 6.4 ± 0.5 Myr. 

Using the above values, we derive My = —14.84 ± 0.38 
and log Lv / Lg = 7.87 ±0.15. The photometric mass is very 
dependent on the form of the IMF. For a Kroupa (2001) 
IMF and a Salpeter IMF (with a lower mass limit of 0.1 
M Q ), we derive masses of 6.6t° 'i x 10 5 and 1.3t° l x 10<5 M o 
respectively. The derived parameters of M82-A1 are given in 
Table H 



Table 4. Summary of derived parameters for cluster M82-A1 and 
its H II region. 



3.4 H 11 region parameters 

We now turn our attention to the properties of the Hn re- 
gion that surrounds the cluster M82-A1. We have already 
established that it is compact with a half-light radius of 
4.5±0.5 pc. From Gaussian fits to the [S 11] emission line pro- 
files, we derive an electron density iV c = 1800^280 cm ~ 3 i as- 
suming a standard electron temperature value of T c = 10 4 K. 
The Hn region surrounding M82-A1 can therefore be classi- 
fied as a high pressure (P/k — 1-2 x 10 7 cm -3 K, considering 



Parameter 



Value 



Cluster 

Radial velocity V T 
F555W 

Half-light radius R c ff 

E(B - V) 

M v 

Log L v /Lq 
Age 

Mass M 

Photon Luminosity 

in Lyman continuum Qo 

H II Region 

Half-light radius R c ff 
Electron density N e 
\ogO/H + 12 
Ionization 
Parameter log U 
Ionized Mass 



+320±20 km s" 1 
17.23±0.10 mag 
3.0±0.5 pc 
1.35±0.15 mag 
-14.84±0.38 mag 
7.87±0.15 
6.4±0.5 Myr 
6.6±g;| x 10 5 Mq 1 
l-3±g;| x 10 6 M Q 2 
7.5±3.0 x 10 50 s" 1 



4.5±0.5 pc 

1800 ±280 cm " 3 

~8.8 
-2.24±0.18 

>5 000 M Q 



1 using a Kroupa (2001) IMF; 2 using a Salpeter IMF. 

the likely range in T c ) compact Hn region. We are also able 
to derive N e for the diffuse gas along the STIS slit in re- 
gion A and find a mean value of ~ 1000 cm -3 , suggesting 
that the ionized gas in region A is at high pressure, at least 
the region probed by the optical [Sll] lines. O'Connell & 
Mangano (1978) derived a similar value of iV = 1800 cm -3 
for regions A and C from the [S 11] ratio. 

Our iV measurements are far higher than those of 
Forster Schreiber et al. (2001) who measured N D in the star- 
burst core using observations obtained with the Short Wave- 
length Spectrograph (SWS) onboard the Infrared Space Ob- 
servatory. The SWS aperture covers an area of 14" x 20" 
and they find that the line ratios of the fine structure lines 
of [S in], [Nem] and [Arm] all lie at the low density limit; 
they adopt N c — 300 cm -3 . The large difference in density 
between the STIS and ISO measurements may be due to 
sampling different environments. We consider this in more 
detail in the next section when we discuss the evolution 
of the M82 starburst as revealed by the study of Forster 
Schreiber et al. (2003). 

We can obtain an approximate oxygen abundance for 
the Hn region using techniques that rely on the relative 
strengths of the strong emission lines when T e is unknown. 
We will use two recent empirical methods. First we follow 
Kewley & Dopita (2002) and use their equation (7) to de- 
rive O/H from the de-reddened [Nil]/ [On] ratio of 2.2. We 
find \ogO/H + 12 = 9.23 based on a solar O abundance 
of logO/H + 12 = 8.9. Thus the O abundance using this 
method is ~ 2x solar. Next, we use the method of Pettini & 
Pagel (2004) which is based on the ratios of [N ll]A6583/Ha 
(N2) and [Om] A5007/H/3 to [N ii]A6583/Hq (OSN2) and 
is insensitive to reddening. Using the relations given by Pet- 
tini & Pagel, we derive 12 + log O/H = 8.76 (jV2 = -0.25) 
and 8.76 (03N2 — —0.087) based on a solar O abundance 
of 8.69. Thus this method gives an O abundance of ~ 1.2 x 
solar. Finally, if we take the derived iV and adopt a value 
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Figure 9. Comparison of the de-reddened (E(B — V) 
50 Myr. The flux is in units of 10 — 14 erg s _1 cm -2 A - 
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0.64) STIS spectrum of M82-F with a synthetic BC03 spectrum for an age of 



of T c = 10 4 K, we derive 12 + log O/H = 8.77. Although 
these methods are approximate, they indicate that the oxy- 
gen abundance of the Hn region surrounding M82-A1 is 
close to solar, and that our assumption that T c = 10 4 K is 
reasonable. 

For comparison, Forster Schreiber et al. (2001) derive 
Ne, Ar and S abundances from their SWS dataset. They find 
that the abundances of Ne and Ar are slightly above solar 
and S is 0.25 solar. In a recent paper, Origila et al. (2004) 
have obtained abundances for red supergiant stars in the 
nuclear region of M82 and find a solar oxygen abundance. 

We can also derive the ionized mass of the Hn region 
from the observed total Ha flux listed in Table [5] Using the 
adopted value for E(B — V), a distance of 3.6 Mpc, T c — 10 4 
K, and -/V e = 1800 cm -3 , we find an ionized mass of > 5 000 
Mq . This value is a lower limit because the close agreement 
between the ages derived from the cluster spectrum and the 
equivalent width of nebular Ha strongly suggests that the 
Hn region is radiation-bounded. 



4 OTHER CLUSTERS 

We fit population synthesis models to the STIS spectra of 
the other four observed clusters using the method described 
above for M82-A1. The quality of the spectra, however, is 
lower than for Al. In Table|KJ we present a summary of the 
derived parameters for the four clusters and compare them 
with previous work. 



4.1 Clusters F & L 

The location and radial velocity of clusters M82-F and L 
(Gallagher & Smith 1999) places them within the 880 pc 
wide gaseous ring surrounding the starburst core of the 
galaxy (Telesco et al. 1991). The low velocity resolution of 



the STIS spectra does not allow us to investigate the kine- 
matics of the clusters. However, the relatively high S/N and 
long wavelength coverage of the spectra allow us to use syn- 
thesis techniques to ascertain the spectral parameters of the 
clusters. 

Figure|5|shows the de-reddened, near-infrared spectrum 
of cluster F compared to the best fit BC03 model. (The 
SB99 models do not offer good spectral resolution here.) The 
quality of the fits is determined by a weighted \ 2 statistic. 
We give more weight to spectral regions containing the Call 
triplet and the broad molecular TiO and VO features. We 
find that clusters F and L have ages of 50I35 Myr and 65I35 
Myr respectively, with E(B - V) values of 0.64 ± 0.10 and 
1.87 ± 0.39. The largest source of uncertainty is the fact 
that the spectral ageing of ~ 20 — 100 Myr old populations 
in the 5000 — 10000A wavelength range can be mimicked by 
a decrease in reddening. 

Cluster F is associated with Ha and [N 11] A6585 line 
emission. Much of this is diffuse and kinematically unrelated 
to the cluster, but we are unable to properly model this 
background because it appears to be highly variable. An 
upper limit to the emission line equivalent width of Ha in 
the spectrum of F is -0.3 ± 0.2A, from which the SB99 
models predict a lower limit to the cluster age of 22 Myr, 
consistent with the estimate above. 

Our age estimate for cluster F compares favourably with 
that of Gallagher & Smith (1999) who find a value of 60 ±20 
Myr (Table 01. However, they derive a reddening value of 
~ 0.9, somewhat in excess of our estimate. This might be 
an effect of patchy extinction in the vicinity of cluster F com- 
bined with the lower spatial resolution of their ground-based 
spectrum. We do notice that the slope of the cluster spec- 
trum becomes progressively redder towards the outskirts of 
the cluster, which might be taken as evidence for a non- 
uniform dust screen or internal extinction gradients within 
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Figure 10. Comparison of the de-reddened (E(B — V) = 0.39) STIS spectrum of M82-B1-1 with a synthetic BC03 spectrum for an age 
of 350 Myr. The flux is in units of 10 — 14 erg s — 1 cm -2 A -1 
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Figure 11. Comparison of the de-reddened (E(B — V) = 0.97) STIS spectrum of M82-B2-1 with a synthetic BC03 spectrum for an age 
of 350 Myr. The flux is in units of 10 — 14 erg s — 1 cm -2 A -1 



the cluster itself. This is consistent with the idea that M82- 
F is seen through a hole in the projected dust distribution 
of the galaxy (Gallagher and Smith 1999). As shown in Ta- 
ble|3 McCrady et al. (2005) find an even lower value for the 
reddening towards F based on ACS/HRC photometry. 

Clusters F and L have comparable ages, which are dis- 
tinct from those in regions A and B. Cluster L has consider- 
ably higher extinction despite being near to F in projection. 
Near infrared (1-2/xm) imaging does not reveal any other 



clusters in their vicinity (e.g. McCrady et al. 2003), so they 
do not seem to be members of a significant intermediate age 
subsystem unless the foreground optical depth in dust is 
very large. In addition, Gallagher & Smith (1999) find that 
clusters F and L have quite different radial velocities of +35 
and +131 kms" 1 respectively. 

Finally, we measured the spatial FWHM of the cluster 
F in different spectral bins to look for any systematic de- 
crease in the cluster size towards redder wavelengths. After 
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Table 5. Derived parameters for clusters M82-F, L, Bl-1 and 
B2-1 and comparison with previous work 



Cluster 


Parameter 


1 his Faper 


Previous Work 


F 


Age (Myr) 




60 ± 20 1 








40-60 2 




E(B - V) 


0.64 ±0.10 


0.9 ± 0.1 3 








0.19-0.37 2 


L 


Age (Myr) 




~ 60 1 




E(B - V) 


1.87 ±0.39 


> 2 3 


Bl-1 


Log (age) (yr) 


8.54 ±0.4 


8.9 4 








8.86 5 




E(B — V) 


0.39 ±0.16 


0.13 4 








0.04 5 


B2-1 


Log (age) (yr) 


8.54 ±0.8 


9.7 4 








9.76 5 




E(B-V) 


0.97 ±0.51 


0.22 4 








0.00 5 



1 Gallagher & Smith (1999); 2 McCrady et al. (2005); 3 Smith & 
Gallagher (2001); 4 de Grijs ct al. (2001); 5 de Grijs et al. (2003). 



correcting for the known wavelength dependent variation of 
the G750L PSF, we find no significant secular changes in the 
size of the cluster. We estimate an average spatial FHWM of 
3.6 pixels or about 183 mas. This compares favorably to the 
average size of the cluster measured from ACS/HRC F814W 
images by McCrady et al. (2005) . 



4.2 Clusters Bl-1 and B2-1 

In Figures 1101 and 1111 we present the spectra and best-fit 
models for the two brightest clusters in regions Bl and B2. 
Both clusters show strong Balmer absorption lines and a 
deep Balmer jump, which are signatures of an intermediate 
age star cluster. The individual cluster spectra are typical of 
the integrated spectrum of region B and most of the rest of 
the main body of M82 outside the starburst core (described 
in de Grijs et al. 2001). (Seen in a distant galaxy, this would 
be called a post-starburst or 'E±A' type spectrum.) From 
our fitting routines, we derive log[age(yr)] = 8.54 ± 0.4 for 
Bl-1 and log[age(yr)]= 8.54±0.8 for B2-1. The significantly 
higher uncertainty for the second cluster is a combination 
of the larger extinction (E(B-V)= 0.97 ± 0.51, as opposed 
to E(B-V)= 0.39 ± 0.16 for Bl-1) as well as a lower average 
S/N. 

These results are generally consistent with those from 
ground-based spectrophotometry and HST imaging of re- 
gion B, which contains a large number of intermediate age 
SSCs (Marcum & O'Connell 1996, de Grijs et al. 2001), 
probably associated with an intense star formation burst 
during the last close passage of M81 and M82 about 1 Gyr 
ago. However, as shown in Table |SJ our ages are younger 
and extinctions higher than those for these clusters based on 
WFPC2 BVI photometry by de Grijs et al. (2001, 2003). We 
believe our spectroscopic ages are more robust, owing to the 
difficulty of doing aperture photometry in such a complex 
region. 



5 DISCUSSION 

The combination of high spatial resolution HST/STIS opti- 
cal spectra and improved spectral synthesis model spectra 
have allowed us to derive numerous parameters for the young 
massive cluster M82-A1 and its Hn region in the starburst 
core of M82. These parameters are summarised in Table [3] 



5.1 M82-A1 and the Evolution of the M82 
Starburst 

We now discuss the age and location of cluster M82-A1 in 
terms of the history of the central starburst. As discussed 
in Section 1, O'Connell et al. (1995) put forward various 
arguments to suggest that region A (with C and E) is the 
part of the starburst core that is least obscured along the 
line of sight. Cluster M82-A1 is located 130 pc north-east 
of the 2/im nucleus, slightly beyond the eastern edge of the 
100 pc radius ionized arc identified by Achtermann & Lacey 
(1995) from [Nell] observations. Its measured radial velocity 
of +320 ± 20 kms" 1 agrees well with the position- velocity 
maps of Achtermann & Lacey (1995). The cluster thus ap- 
pears to be spatially and kinematically associated with the 
ionized gas seen in the starburst core of M82. 

The evolutionary synthesis models of Forster Schreiber 
et al. (2003) reproduce the observed characteristics of the 
M82 starburst if there were two short duration bursts occur- 
ring 8-15 and 4-6 Myr ago. The first burst was concentrated 
in the central few tens of parsecs and the second episode oc- 
curred mainly in the circumnuclear regions. On the basis 
of its position and age of 6.4 Myr, we associate M82-A1 
with the more recent burst. Melo et al. (2005) catalogue 86 
clusters in region A and find this region has the highest pro- 
jected surface density of young clusters. The morphology of 
region A (cf. Fig. suggests the clusters are coeval. This is 
strengthened by the fact that the bright base (in projection) 
of the superwind is centred on regions A and C (O'Connell 
& Mangano 1978; Shopbell & Bland-Hawthorn 1998; Gal- 
lagher et al., in prep.), indicating that the clusters must be 
of a similar age since the cluster winds are very time depen- 
dent. If this is the case, then region A must correspond to a 
particularly intense episode of star formation ~ 6 Myr ago. 



5.2 The M82-A1 H II Region 

We now discuss the compact Hn region associated with 
M82-A1. Forster Schreiber et al. (2001) have studied the 
physical conditions of the ISM in the starburst core of M82 
from ground and space-based IR observations. They find 
that a random distribution of closely-packed ionizing clus- 
ters and small gas clouds, separated on average by a few 
parsecs, represents well the starburst regions of M82. They 
derive a constant ionization parameter log U of —2.3 on spa- 
tial scales of a few tens of parsecs to 500 pc, and suggest 
that this uniform degree of ionization indicates a similar 
star formation efficiency and evolutionary stage for the M82 
starburst regions. 

The Hn region around M82-A1 represents a photoion- 
ized nebula on a much smaller physical scale than previous 
observations of Forster Schreiber et al. (2001) which probe 
down to structures typically on 25 pc scales. For comparison 
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with their model, we derive the ionization parameter for the 
M82-A1 Hn region as given by: 



U 



Qo 



4ttR 2 N c 



(1) 



where Qo represents the photon luminosity in the Lyman 
continuum, R is the nebula radius, and -/V e is the elec- 
tron density. Taking the values given in Table |1J we ob- 
tain log U — —2.24 ± 0.18, in excellent agreement with the 
value of —2.3 derived by Forster Schreiber et al. (2001), and 
thus extending the uniformity of the M82 starburst down 
to scales of a few parsecs or individual cluster Hn regions. 
This suggests that the whole M82 starburst region is self- 
adjusting in the sense that the number of ionizing photons 
per atom is constant. 

The value of U derived by Forster Schreiber et al. (2001) 
assumed an electron density of 300 cm -3 , which is far lower 
than the value we find for the M82-A1 Hn region and the 
surrounding medium. We suspect these differences are due 
to sampling different environments; the IR observations of 
Forster Schreiber et al. (2001) probably probe deeper into 
the starburst region, and into a region where the gas is less 
compressed. 

Turning now to the size of the Hn region, we find that 
it is compact with a radius only 1.5 times the cluster radius. 
It has a high pressure of P/k = 1-2 x 10 7 cm -3 K. The ambi- 
ent thermal pressure of the diffuse ionized gas in region A is 
also high with P/k ~ 10 7 cm -3 K. The total ambient pres- 
sure is probably even greater than this. The STIS spectra 
of the diffuse gas in region A show that the Ha emission is 
broad with an average velocity dispersion a v of 45 kms -1 , 
giving Ptuib/k ~ pa^/k ~ 3 x 10 8 cm -3 K. This suggests 
that the turbulent pressure component dominates. It pre- 
sumably originates from the nonthermal energy input from 
supernovae events. The M82-A1 Hn region is therefore in 
an unusually high pressure environment which will control 
its evolution. 

Observations of radio supernova remnants (Muxlow et 
al. 1994; Pedlar et al. 1999) in M82 show that they are com- 
pact (with radii < 4 pc) for their inferred ages; these authors 
suggest that their growth has been inhibited by the dense 
starburst environment in M82. Chevalier & Fransson (2001) 
propose that these sources are evolved supernova remnants 
interacting with a dense interstellar medium; they may have 
recently entered the radiative phase if the surrounding den- 
sity is ~ 1000 cm -3 . We find that the density in region A is 
of this order, and that the compact Hn region surrounding 
M82-A1 has a similar size to the upper limit on the sizes of 
the radio supernova remnants. We thus postulate that the 
compactness of both the SNRs and the M82-A1 Hn region 
is the result of the high ambient pressure in the starburst 
core. 

The continuous injection of energy by winds into the 
interstellar medium results in the formation of an adiabatic, 
pressure-driven bubble (e.g. Weaver et al. 1977; Chevalier 
& Clegg 1985). For a star cluster, the source of energy is 
first stellar winds and then supernovae. For M82-A1, we es- 
timate that its cluster wind has a mechanical luminosity 
imoch ~ 2.5 x 10 40 ergs -1 from the Starburst99 models for 
a Kroupa IMF and the derived mass of M82-A1. We can 
estimate the radius of the shell formed by the cluster wind 



using the standard equation for pressure-driven superbub- 
bles from Weaver et al. (1977): 
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where R is the shell radius, po is the ambient mass density, 
and t is the age. Using the value of L mcc h given above, an 
ambient number density of 1000 cm -3 (Sect. ET41 . and an 
age of 6.4 Myr, the predicted radius is 150 or 100 pc, if 
we assume a thermalisation efficiency of 100% or 10% re- 
spectively. These values are obviously far in excess of the 
measured radius of 4.5 pc, and suggest that the Hn region 
has not developed according to the standard model. 

Oey & Gara'a-Segura (2004) and Dopita et al. (2005) 
consider the effect of the ambient interstellar pressure on the 
dynamical evolution of pressure-driven superbubbles. A high 
ambient pressure will inhibit the growth of the H II region 
and it will stall when the internal pressure of the superbub- 
ble equals the ambient pressure. The internal pressure P is 
given by Weaver et al. (1977) as: 
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This leads to an internal pressure P/k of ~ 3 x 10 9 or 
~ 6 x 10 8 cm -3 K, assuming thermal efficiencies of 100 and 
10% respectively. The latter value is reasonably close to our 
estimate of the total (thermal and turbulent) ambient inter- 
stellar pressure of ~ 3 x 10 8 cm" 3 K. While this is a rough 
estimate, and does not take into account the fact that the 
bubble will not abruptly stop expanding when it stalls, it 
illustrates that the high ambient pressure is probably suffi- 
cient to modify the evolution of the H II region surrounding 
the cluster M82-A1. 

Using Eqn. (3), Dopita et al. (2005) show that all stalled 
Hll regions will have a common ratio of L moc h/ p R? , and 
since I/ me ch scales roughly with the flux of ionizing photons 
Qo, then all stalled H II regions should have the same ioniza- 
tion parameter as given by Eqn. (1). This relationship may 
help to explain the uniformity of the ionization parameter 
on many scales in the starburst regions of M82 discussed 
above and by Forster Schreiber et al. (2001). 

Another important possibility to consider is whether the 
cluster wind from M82-A1 has developed according to the 
standard adiabatically expanding wind model. Canto, Raga 
& Rodriguez (2000) present an analytic mass-loaded cluster 
wind model for a dense cluster of massive stars. They find 
that there are two solutions depending on the properties of 
the surrounding ISM; the standard supersonic solution or 
a subsonic solution if the environment is at high pressure. 
The latter solution may apply to M82-A1. The supersonic 
wind solution can also be inhibited if the cluster wind is suf- 
ficiently dense to promote radiative cooling (Silich, Tenorio- 
Tagle & Munon-Tunon 2003). This occurs if the cluster is 
massive but compact and a supernebula can then be formed 
(Tenorio-Tagle et al. 2005). From their figures, we find that 
the cluster M82-A1 is not in the regime where catastrophic 
cooling will occur. 

Ultimately, it is hard to understand how the observed 
M82 galactic superwind is generated if the environment sti- 
fles its development because of the high ambient pressure. 
Chevalier & Fransson (2001) consider this problem in rela- 



16 L. J. Smith et al. 



tion to their prediction that the compact radio supernovae 
in M82 are radiative and thus are unable to provide suffi- 
cient energy to drive a galactic wind. They hypothesise that 
galactic winds are driven only by supernovae occurring in a 
low density, hot medium. Another way around this problem 
is to conjecture that there is a pressure gradient in the cen- 
tral region of M82 such that the outer clusters are in a lower 
pressure region. The winds from these clusters will then be 
able to develop and may produce highly-collimated outflows, 
as envisaged in the interacting wind model of Tenorio-Tagle, 
Silich & Munon-Tunon (2003). 

Whatever the precise reason for the compact nature of 
the Hn region surrounding M82-A1, it is clear that the un- 
usually high ambient pressure in the M82 starburst environ- 
ment must play an important role. 



6 SUMMARY AND CONCLUSIONS 

We have presented optical HST/STIS spectroscopy of five 
massive clusters in M82. Our best data are for cluster M82- 
Al and its environment in region A of the starburst core. 
This region contains a large complex of young clusters which 
appear to be responsible for a major component of the su- 
perwind. We have sought to derive the properties of M82-A1, 
and thereby clarify the age of the recent starburst, and to 
investigate the nature of its ionized environment. 

The STIS G430L spectrum of M82-A1 is dominated by 
Balmer absorption lines and a weak Balmer jump suggesting 
a young age. We measure a mean heliocentric velocity of 
+320 ±20 kms -1 . This value and its position 130 pc north- 
east of the 2/im nucleus suggests that M82-A1 is spatially 
and kinematically associated with the starburst core. 

Using the evolutionary synthesis code Starburst99 (Lei- 
therer et al. 1999) which includes a high resolution spectral 
library to derive both the age and reddening of M82-A1, we 
find that we can only match the spectral energy distribu- 
tion and depth of the Balmer jump for a narrow age range 
of 5-7 Myr and a reddening of E(B - V) = 1.5-1.2 mag. 
The independent BC03 models give similar results. 

The STIS G750M spectrum shows that M82-A1 is sur- 
rounded by a compact Hn region. We have recovered the 
intrinsic nebular spectrum by subtracting the best-fitting 
synthetic spectrum and find that the nebular reddening of 
E(B — V) = 1.38 ± 0.15 mag agrees well with the value de- 
termined from the cluster SED. We adopt a mean of the two 
methods to give E(B ~V) = 1.35 ± 0.15 mag. The observed 
equivalent width of the Ha emission line provides an upper 
limit to the age of 6.7 Myr, in agreement with the depth 
of the Balmer jump. We adopt an age of 6.4 ± 0.5 Myr for 
M82-A1. 

Comparison with the models of Forster Schreiber et al. 
(2003) suggest that the formation of M82-A1 is associated 
with the most recent burst of star formation which took 
place «4-6 Myr ago. The overall morphology of region A and 
the fact that a major component of the superwind is centred 
on this region suggests that the clusters are coeval. If this 
is the case, then region A must correspond to a particularly 
intense episode of star formation. 

We have used archive HST images to derive the radius 
and mass of M82-A1 via photometry. We find that the clus- 
ter is elliptical with a mean effective radius of 3.0 ± 0.5 pc 



and a minor/major axis ratio of 0.72 ± 0.05. The cluster 
effective radius appears to decrease towards shorter wave- 
lengths. We interpret this behaviour as being due to dust 
in the surrounding H II region rather than mass segregation 
because it seems difficult to produce a radial colour gradi- 
ent in a young cluster. We derive a mass of 0.7-1. 3xl0 6 M 
for M82-A1 depending on whether we adopt a Kroupa or 
Salpeter IMF. 

We have sought to understand the nature of the H n re- 
gion surrounding M82-A1. From HST/PC images, we have 
derived an effective radius of only 4.5 ±0.5 pc and find that 
it is highly elliptical with a minor/major axis ratio of 0.57. It 
is also of high density; the [S n] doublet in the STIS G750M 
spectrum gives an electron density of 1800 ± 300 cm~ or a 
thermal pressure of P/k = 1-2 x 10 7 cm~ 3 K, for T c =0.5- 
lxlO 4 . The diffuse ionized gas along the STIS slit has a 
slightly lower density which, together with a measured high 
turbulent velocity dispersion of 45 kms -1 , suggests that 
M82-A1 is immersed in an unusually high pressure environ- 
ment. We also find that the oxygen abundance of the Hn 
region is close to, or slightly higher than, solar. 

We find that the ionization parameter log U of the H II 
region is —2.24 ±0.18 in excellent agreement with the value 
of —2.3 derived by Forster Schreiber et al. (2001) for scales of 
a few tens to 500 pc. Our measurement thus extends the uni- 
formity of the M82 starburst down to the scale of individual 
cluster Hn regions, and suggests that the entire starburst 
is self-adjusting with a constant number of ionizing photons 
per atom. 

We have sought to understand the compactness of the 
Hn region around M82-A1, given that it is much smaller 
than the size predicted for its age according to the standard 
model for a pressure-driven superbubble. We find that the 
high pressure of the environment may have caused the bub- 
ble to stall or follow a different wind solution to the standard 
supersonic adiabatic one. We suggest that models should be 
developed which account for the evolution of a cluster wind 
in a high pressure environment. 

We also derive spectroscopic ages and reddening for four 
other clusters in M82 (Bl-1, B2-1, F and L). These are con- 
sistent with earlier studies and demonstrate that star for- 
mation activity, sufficiently intense to produce super star 
clusters, has been going on in M82 during the past Gyr, 
perhaps in discrete and localized episodes. 
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